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Abstract. We have developed a numerical model for the temporal evolution of particle and photon 
spectra resulting from nonthermal processes at the shock fronts formed in merging clusters of 
galaxies. Fermi acceleration is approximated by injecting power-law distributions of particles during 
a merger event, subject to constraints on maximum particle energies. We consider synchrotron, 
bremsstrahlung, Compton, and Coulomb processes for the electrons, nuclear, photomeson, and 
Coulomb processes for the protons, and knock-on electron production. Broadband radio through y- 
ray light curves radiated by nonthermal protons and primary and secondary electrons are calculated 
both during and after the merger event. Using ROSAT observations to establish typical parameters 
for the matter density profile of clusters of galaxies, we find that merger shocks are weak and 
accelerate particles with relatively soft spectra. Our results suggest that only a minor contribution to 
the diffuse extragalactic y-ray background can originate from cluster merger shocks. 



INTRODUCTION 

In the hierarchical merging cluster scenario, clusters grow in mass by accreting nearby 
clusters. Approximately 30-40% of galaxy clusters show evidence of substructure in 
both the optical [1] and X-ray wavelengths [2]. Velocity differences between the ob- 
served structures is ~ 1000-2000 km s _1 . With gravitational forces driving the inter- 
action between the two systems, cluster mergers are consistent with highly-parabolic 
orbits. Typical sound speeds within the intracluster medium (ICM) are ^1000 km s~ ! , 
so shocks will form at the interaction boundary of the two systems. 

Shock fronts that form in the ICM as a result of a cluster merger event are thought to be 
associated with the cluster radio relics, which are diffuse emission found on the cluster 
periphery with no known optical counterpart. The shock compression will orient any 
existing cluster magnetic field into the plane of the shock. Radio relics are characterized 
by highly organized magnetic fields with field strengths in the ~1 /iG range with linearly 
polarized field lines in the vicinity of the shock [3]. The shock front will accelerate a 
fraction of the thermal particles within the ICM by first-order Fermi acceleration. 

It has been proposed that cluster mergers are the dominant contributor to the diffuse 
y-ray background [4]. Some unidentified EGRET sources are claimed to be associated 
with y-ray emission from galaxy clusters [5]. Excess EUV emission from Coma can be 
explained by nonthermal electrons accelerated at merger shocks [6]. Variations in radio 
surface brightness will result from superposition of cluster emissions [7]. 



MODELS 



We present the results of a computer code designed to calculate the time-dependent par- 
ticle distribution functions evolving through adiabatic and radiative losses for electrons 
and protons accelerated by the first-order Fermi process at the cluster merger shock [8]. 
The model calculates the shock speed from the shock formed in a cluster merger event 
by the interaction speeds expected from two point masses interacting under their mutual 
gravity. The point masses are assumed to be on elliptical orbits whose onset is initiated 
by the collapse and merger of primordial density fluctuations. The collision is a result 
of the two bodies "falling out" of the Hubble flow and onto a nearby cluster. Expected 
total orbital energies are E tot ~ 10 63 -10 64 ergs. Expected collision velocities between a 
1O 15 M and a 1O 14 M mass cluster range from ~1 800-3500 km s 1 . The electron and 
proton distribution functions originate from a momentum power-law injection spectrum. 
In terms of particle kinetic energy E e _ p = (y e , p — l)m e p c 2 , the injection function is 

Qe, P ( E e,p,t) = Q% p c s{t) ~ l [ E e, P (E e:P + 2m e , p c 2 )Y^ (E e , P + m e , p c 2 ) e~T&k , (1) 

where an exponential cutoff E m&x has been applied and is determined by the maximum 
energy associated with the available time since the beginning of the merger event, by a 
comparison of the Larmor radius with the size scale of the system, and by a comparison 
of the energy-gain rate through first-order Fermi acceleration with the energy-loss rate 
due to adiabatic, synchrotron and Compton processes. The injected particle function is 
normalized through the normalization constant Q Q e which is determined by 
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dE E e , p Q e , p (E,t) = ^Asn^mpv] . (2) 



where r\ e , p is an efficiency factor, v s is the shock speed A s is the area of the shock 
front; E mSLX and E min are the maximum and minimum particle energies; is an 
enhancement factor to account for the presence of ions heavier than Hydrogen. We 
assume an efficiency factor r\ e p = 5% for both protons and electrons. With this method 
the total energy deposited into nonthermal particle production is r\ e ^ p E toV The particle 
density (n) is calculated according the /3 model, and averaged over the shock front. 

The time evolving particle spectrum is determined by solving the Fokker-Planck 
equation in energy space for a spatially homogeneous ICM, given by 

^4|, [I , (£ , 0W ( £ , 0] _^ [fi , ol(£ , (M£ , 0]+e ( £ , _j^^M. 

(3) 

The quantity E tot (E,t) represents the total synchrotron, Compton, Coulomb, and adia- 
batic energy-loss rate for electrons, and the sum of the Coulomb and adiabatic energy- 
loss rates for protons. Both protons and electrons are subject to diffusion in energy space 
by Coulomb interactions. The protons experience catastrophic losses due to proton- 
proton collisions on the timescale x n , proton-y collisions on the timescale x p y, and spatial 
diffusion from the cluster on timescale % d . Spectra of secondary electrons and positrons 
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FIGURE 1. Light curves at various observing frequencies produced by a shock formed in a merger 
between 1O 14 M and 10 15 M p clusters that begins at z ( =0.3 and is evolved to the present epoch (t=3.42 
Gyr). All light curves are for a magnetic field strength of B=l .0 jiG except the 15 MHz light curve is also 
calculated with a magnetic field strength of B=0.1 fxG (dotted curve). Radio light curves in Jansky units 
are given at 15 MHz, 30 MHz, 74 MHz, 110 MHz, and 1.4 GHz on the left panel, and light curves in 
energy flux units are given at 40-80 keV, > 100 MeV, > 1 GeV and > 100 GeV in the right panel. 



are calculated from pion-decay products, and are subject to the same physical processes 
as the primary electrons. Q(E,t) is the particle injection function. 

The synchrotron, Compton, bremsstrahlung, and pion-decay y-ray spectral compo- 
nents are calculated from the particle spectra following the methods described by [8]. 
We use a standard parameter set with a mean ICM number density n Q = 10~ 3 cm~ 3 and 
a uniform cluster magnetic field of B = 1 .0 /lG. Thus, V (t) =A(t)v s t. 



RESULTS 

Fig. 1 shows the light curves at various energies for a cluster merger shock between 
a 10 14 M Q and 10 15 M Q with radii of 1.5 Mpc and 0.75 Mpc, respectively, with a 
central gas number density of n =10~ 3 , and a magnetic field of B=l /iG is assumed. 
The assumed core radius for the dominant cluster profile is 250 kpc, and /3=0.75. We 
assume a cosmology defined by H = 70 km s" 1 Mpc" 1 , and (Q ,Q^,Q A )=(0.3,0.0,0.7). 
We assume the shock forms at a redshift of z—0.3. The light curves exhibit a similar 
characteristic independent of frequency with the peak luminosity occurring when the 
centers of mass of the two clusters pass. The emission slowly decay and approach a 
plateau at times t > t acc when particle injection has stopped. The rate of decay slows with 
decreasing energy with the lowest energies exhibiting the slowest decay rates. The late 
time plateaus are from n° y-rays for photons with E>1 GeV, and secondary electrons 
for the radio energies. With these expected fluxes, it is unlikely that the more than a 
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FIGURE 2. Calculations of the hardest particle injection spectral indices s min formed in cluster merger 
shocks as a function of the dominant cluster mass, for various values of r c and j3 . The mass of the merging 
cluster is assumed to be 10 14 M Q . Our adopted values are plotted by the solid line, and dot-dashed line 
shows s min plotted against the dominant cluster mass for the average values of r c and j3 obtained by 
ROSAT observations of 45 Abell clusters [9]. 



few of the isotropic unidentified EGRET sources can be attributed to the radiation of 
nonthermal particles produced in cluster merger shocks. 

The diffuse extragalactic y-ray background is a featureless power law with photon 
index of 2.10(±0.03). Fig. 2 shows the hardest particle injection spectral index, s min , 
for a cluster merger shock as function of the mass of the dominant cluster with varying 
values of core radii, r c , and /3's. The particle spectral indices expected from the average 
values of r c and /3 obtained from ROSAT observations of 45 Abell clusters [9] are 2.2- 
2.4. Only the dark matter profiles with strong central peaks produce particle indices 
<2.1. This discrepancy in the particle indices suggests that nonthermal radiation from 
cluster merger shocks can make only a minor contribution to the diffuse extragalactic 
y-ray background unless dark matter halos contain strong central peaks. 
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